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Abstract 

Cosmological models where spatial sections are the Poincare dodeca- 
hedral space D have been recently invoked to give an account of the lower 
modes of the angular anisotropics of the cosmic microwave background. 
Further explorations of this possibility require the knowledge of the eigen- 
modes of the Laplacian of D. Only the first modes have been calculated 
numerically. Here we give an explicit form for these modes up to arbi- 
trary order, in term of eigenvectors of a small rank matrix, very easy to 
calculate numerically. As an illustration we give the first modes, up to the 
eigenvalue — k(k + 2) for k — 62. These results are obtained by application 
of a more general method (presented in a previous work) which allows to 
express the properties of any eigenfunction of the Laplacian of the three 
sphere under an arbitrary rotation of SO (4). 

1 Introduction 

There is a long time interest for cosmological models where space is multi- 
rather than simply-connected (see a review in Recently, |10| claimed 
that a peculiar model may give an account of the two first moments in the 
angular power spectrum of the anisotropics of the cosmic microwave back- 
ground. This model involves the Poincare dodecahedral space D, whose 
universal covering is the three sphere S 3 . The calculations involved the 
first eigenmodes [of the Laplacian] in D, which were estimated by numer- 
ical methods. However, to make predictions beyond the first moments, 
and to check non diagonal terms in the correlation matrix, the knowledge 
of a greater number of modes is required. 

The eigenmodes of a spherical space S 3 /T are the an eigenmodes of 
S 3 which are conserved by all the rotations of P. In general, they remain 
unknown. Their number, i.e., the multiplicity of the eigenspaces, was cal- 
culated by |1] as a function of the eigenvalue of the Laplacian. Recently, 
[HI (see also [S]) have provided analytic calculations of these modes for 
Lens space and Prism space. The same results were obtained later with 
a different method by [jj. Here we provide a new and efficient method 
to calculate the modes of the dodecahedral space. In 0, we constructed 
a special basis (B3) for the modes of S 3 , which allowed to calculate ex- 
plicitely their behaviour under any rotation of SO (4). These results are 
recalled in section 2. They are then applied to D (section 3), which leads 
to the derivation of its eigenmodes. 
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2 Rotation properties of the eigenmodes 

of S 3 

The eigenvalues of the Laplacian A of S 3 are of the form \ k = — k (fc + 2), 
where k £ IN + . For a given value of k, they span the eigenspace V k of 
dimension (k + 1) 2 . The eigenmodes of the dodecahedral space D = S 3 /G 
are those eigenmodes of S s which remain invariant under all holonomy 
transformations of G. To be so, it is necessary and sufficient that they 
are invariant under the two generators of G, that we call g±. 

We will use a specific basis B2 of V fc , whose properties were detailed 
in This basis was originally introduced by p^. A real version of it has 
been used by |H] to find the eigenmodes of lens and prism spaces. It is 
generated by a set of (k + l) 2 functions: 

B2=(T k;mi , m2 ), mi, m 2 = -k/2...k/2, (1) 

where mi and 7712 vary independently by entire increments. We know ([!]) 
that D has eigenmodes only for fe even, that we assume from now. This 
ensures that mi and 7772 are entire. Hereafter, all sums involving mi or 
r7i2 will cover the range — fe/2..fe/2. 

These functions are naturally adapted to toroidal coordinates to de- 
scribe S s : {x,9,4>) spanning the range < x < 7r /2, < 8 < 2ty and 
< (j> < 2n. They are conveniently defined (see 8 for a more complete 
description) from an isometric embedding of S 3 in IR 4 (as the hypersurface 
xeM 4 ; I x |= 1): 

x° = r cos x cos 
x 1 = r sin x cos <f> 
x 2 = r sin x sm <t> 
x 3 = r cos x sin # 

where (a; M ), p = 0, 1,2,3 is a point of IR 4 . 

In we gave the explicit expression of B2: 

T k;mi , m2 {X) = C fe;mi , m2 [cosx e B ] 1 [sinx e l0 ] m P^ m ' e) [cos(2x)] ■ (2) 
We wrote £ = mi + 7712, m = 7772 — mi for simplicity; p^* 1 ' 1 ' is a Jacobi 
polynomial. Normalization leads to C k[mi ,m 2 = \J fJ/\Xl% jt/'-mj! ' 

Note that the basis functions Tk :mi ,m 2 have also been introduced in [5] 
(p. 253), with their expression in Jacobi Polynomials. They are the com- 
plex counterparts of those proposed by ||] (their equ. 19). The variation 
range of the indices mi , 7772 here is equivalent to their condition 

\£\ + \m\<k, £ + m = k, mod 2, (3) 
through the correspondence £ — mi + 7772, m = 7772 — mi. 

2.1 Complex null vectors and roots of unity 

Our calculations require the use of an other basis of V k which was intro- 
duced in |S] (and apparently ignored before in the literature): 

S3 = ($fj), I,J = 0..k. (4) 

This basis may only be used when k even, which is assumed here. Here- 
after, all sums involving i,j, I or J will cover the range 0, k. 
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The basis functions are defined as <J?f j : $jj(X) = (X ■ Nu) k , from 
the (fc + l) 2 null vectors of C 4 , 

Nu = N(Ia, Jet) — (cos la, i sin Ja, i cos Ja, sin/a). (5) 

Here the dot product extends the Euclidean [scalar] dot product of 1R 4 
to its complexification C 4 . (Complex) null vectors are defined as having 
zero norm. The angle a = -j^j is the argument of the (fc + l) th complex 
roots of unity. Those are the powers p 1 of 



p=e k + 1 = cos a + i sin a. (6) 

We proved in [B] that B3 form a basis of V k (when fc is even) and we 
gave the transformation formulae with B2. First we defined Tk-m 1 ,m 2 = Pk;m 1 ,m 2 Tk- mi ,m 2 , 
with 

p = ^1 

k; mi ,m 2 (fc/2 _ mi) , (fc/2 + mi) , (fc + 1)2 Ck ^ i m2 

2- k 7T fc! (fc + 1)~ 5/2 

~ + m 2 )! (fc/2 - m 2 )!(fc/2 + mi)! (fc/2 - mi)! ' 

Then the transformation reads 

1 fc 

-7- _ ____ \^ I(m 1 +m 2 )-J(m 2 -m 1 ) 

2k, mi ,m 2 - 7T 1)2 2^ P ®U' ( -'> 

V ; I,J=0 

fc/2 

ifc _ \ ,7- -I(mi+m 2 )+J(ni2-mi) 

yjj — -lk:m 1 ,m 2 p 

vai ,m 2 =—k /I 



2.2 Rotations in H 4 

The main interest of the basis B3 lies in the possibility to calculate ex- 
plicitly its rotation properties under an arbitrary rotation g GSO(4). The 
action of a rotation g £ SO(4) on a function is noted R 9 : R 9 /(x) = f(gx). 
Applied to the basis functions, it may be developed as 

fc 

R 9 : $ 7J i > R 9 $ 7J : = ^ 1J { gx ) = G$M (8) 

ij=0 

The coefficients were calculated in as 

G " = 7]rXT^ E P- I(A+S " fe) P- J(A - S) W A V B , (9) 

^ + A,S=0 

where we defined W = (-^r), V = (^), A =< Ql a Qr ■ nu >, 
A' =< Ql a Q R - nu >, B =< Ql P Qr ■ nu >, V =< Ql S Qr ■ nu > 
. Note that AA' = BD. 

In these formulae, the two unit quaternions Ql and Qr represent the 
rotation g. The four basis quaternions are written j^, j — 0, 1, 2, 3, with 
j° — 1. A quaternions is developed as q = q M j^. For a complex quater- 
nions, the q M are complex numbers. Terms like Ql a Qr simply denote 
the quaternionic product. The brackets indicate the quaternionic scalar 
product, the bar the quaternionic conjugate. The null complex quater- 
nions nu represent the null vectors Nu introduced above (see details in 
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0). Also, we introduced the following peculiar null complex quaternions 
a s 1 + i js, P = ji — i j% = (1 — i js) j\ and 5 = —ji — i fa. They have 
zero norm and obey the properties 

<a-n IJ >=p I , <a-n IJ >=p~ I , < (3-riu >= p J , <S-nij>=p~ J . 

The coefficients G/'j completely define the transformation properties 
of the basis functions under any element g of SO(4). 



3 Dodecahedral space 
3.1 Generators 

The Poincare Dodecahedral space has two generators g±, acting asm 
g±x, for any point of S 3 represented by the unit vector x of the embedding 
space ]R 4 . In a certain basis of IR 4 , they are expressed (Weeks, prived com- 



munication) as the two matrices g± 
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±1/2 





c 


-C 


V o 
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with C ; 
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It follows their (left action) complex matrix forms as G± 
whose action is defined by 



' w 


iZ ' 


H-> G± 


' w 


iZ ' 


iZ 


w 


iZ 


w 



TT 7" 0,-3 

W = x ± ix , 



c±iC 
±1/2 



ia; G C. 



±1/2 
c- iC 



These two operators have the same eigenvalues, namely A = e™/ 5 and 
A* (star means complex conjugation). 

The corresponding (unit) quaternions Q± act by left action only (as 
<lx i— > Q± Qx, where q x is the unit quaternion associted to the point x of 
S 3 ). They are given by Q± = c ± C j 3 ± ji/2. We note A = e j37r/5 the 
quaternionic analog of the eigenvalue A. 

Diagonalisation 

Calculations are much easier if we adopt a different basis where one 
generator, say g+, transform the complex coordinates W = x° ± ix 3 and 
Z = x 1 ± ix 2 , of an arbitrary vector a; € IR 3 , by scalar (although com- 
plex) multiplication. To find such a basis is equivalent to diagonalize the 
complex matrix G+ , what we have done. We express the results in quater- 
nionic notation. We note u the quaternion which expresses the change of 
coordinates: by definition, A = u Q+ w _1 . Calculations give 



u = 1/2 +j 2 U/C 2 + l/4-C). 

The second generator takes the form 

R^=uQ^ u' 1 = cos(7r/5)-j 3 sin(7r/5) /VE + ji 2 sin(7r/5) /VE. (10) 

Now we continue the calculations with this basis: the two generators are 
A (diagonal) and i?_ (formula 11011 . 



3.2 Invar iance 

The diagonal character of the first generator makes the calculations easy: 
the basis functions Tk;M x ,Mz °f B2 appear to be its eigenfunctions. This is 
analog to the case of a lens space, as seen in 0. Thus, the functions of V k 
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which are g+ invariant are all combinations of the Tk ]mi ,m 2 which verify 
the condition 2m2 = mod fO. This implies fc even (as was derived by 
[U), so that mi and m.2 are entire and 

m2 = mod 5, (ff) 

where we underline as ni2 a value of 7712 verifying this condition. It remains 
to express the invariance condition with respect to the second generator. 
Second generator 

We assume fc even, and we adopt the basis B3 of V . Simple calcula- 
tions lead to estimate 

Au =< (R-)a-nu >= F p'+E p J , A'u =< (R-)a-nu >= F* p^-Ep'J, 



Bu =< {R-) P P ■ nu >= F* p J - E P p 1 , (12) 

with F = cos(7r/5) + i sin(7r/5)/V5, E = 2 sin(7r/5)/- v / 5- 

We insert the relation U = -jj = p I+J in Taking into account the 
properties of the roots of unity, we obtain 

/ a) '\k —i -1 

G%{R-) = l 7 ' p 8 Dlrac V 



(fc + 1) V pI + J-2i - 1' 

where 5 Dlrac holds for S Dirac [(I + J) - (i + j), mod (k + 1)]. 

This expression allows to return to the rotation properties of the ba- 
sis B2 (in fact Th-M 1 ,M 2 rather than Tt-M 1 ,M 2 ), that we express by the 
development 

R.g_T k ; Ml ,M 2 = ^2 r MiMj(9) 2fc;mi,m 2 - (13) 

The coefficients take the form 



r m 1 m 2 \ "» \ " I(M 1 +M 2 ) + J(M 1 -M 2 )-i(m 1 +m 2 )-j (m 1 -m 2 ) <fjj 
L M 1 M 2 - P (k + 1) 2 

a i. j K ' 

(14) 



+(I+J) (-mi+ma+MQ [V - jj J P 



P (fc + 1) 3 (V/+ J -'- ()•) 

Direct calculations allow to evaluate the last sum as 

leading finally to 

Vk/2 — mo ( A 1 \ fe 

| .... X , : - l-^ ) (I -\ 

Ml M 2 P (fc+1) 2 ' ( ' 

Still a new basis 

This suggests us the introduction of still a new basis for V k (only for 
intermediary calculations) : 

f aM2 = Y, p a Ml T k; Ml ,M 2 , « = O.-fc, M 2 = -fc/2..fc/2, (16) 

with inverse formula Tk-,M X ,M 2 = jjij J2 a P°" Ml T a M 2 - Its rotation prop- 
erties are deduced as: 

a (k/2-M 2 ) _ 

R 9 : f aM2 i ► P — J2 P 21 M2 ( A ')" ^ E V_m2 ( 17 ) 
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where the quantities A' and V have to be evaluated for the value J = a — I 
mod (k + 1). 

A first important result appears: the value of a, for this basis, is 
preserved by R g . Using the formula just above, it is straightforward to 
check that this implies, similarly, that R 9 preserves mi in the basis B2. 
This result will allow considerable simplification. To continue, we report 
in the formula above the values of V = A/B and A' (eau WIl . This gives 

R 9 : f aM2 i-> 2J tS 2 T am2 , 

m 2 

with 

(/-a/2) (2M 2 +fc) 

'Cm =Y ~ i ; [F + E p a - 2I ] k/2 - m2 [F* p a ~ 21 -E] k/2+m \ 

2 *■ — ' k + 1 

(18) 

where k/2 — m,2 and k/2 + rri2 take their value between and k. 



3.3 Splitting the eigenspace 

From this formula, it appears the second important result that 

_ I (2M 2 + k) 

iTm 2 = iZ_ = E k + 1 l F + E p- 2I ] k/2 ~ m2 [F * P - 21 - E] k / 2+m * 

(19) 

does not depend on a. This allows to write the rotation formula 

m 2 

It is easily checked that this absence of dependence on a holds also for 
the 7fc :mii ,„ 2 . Finally, we are led to the transformation rule for the basis 
B2: 

R. s _ : Tk;M u M 2 ^2lMl Tk;M 1 ,m 2 , (20) 
in 2 

with 7S2 s iven b y CHI- 

Summary 

We summarize the results obtained, expressed in the basis (7fc ;mi , m2 ): 

• The g+ invariant functions of V k are all combinations of the 7fc ;mi , m2 
which verify the condition m,2 = mod 5. 

• The second generator (g~) preserves the value of mi. This allows 
to consider V k as the direct sum of k + 1 sub vector-spaces V h ' mi 
of dimension k + 1, each preserved by g-. The search for invariant 
functions can thus be made independently for each V k,mi . 



• The rotation coefficients 1191 in V k ' mi do not depend on mi: they 
are identical in each of the V k ' mi . It is thus sufficient to search 
invariant functions in one of them (say, for mi = 0). The other 
functions are given by the k more copies obtained by replacing the 
value mi = by any of the k remaining values. It results that the 
dimension of the space of eigen functions of the dodecahedral space 
is an entire multiple of k + 1, as was already indicated by 



G 



3.4 Eigenvectors 

The search of the modes of the dodecahedral space is reduced (for each 
even value of k) to that of the invariant eigenfunctions of S 3 in y k < m i= l \ 
Invariance with respect to g+ implies that such a function may decom- 
posed in the restricted basis formed by the 7m 1= o,m 2 (we recall that un- 
derlining means M 2 = mod 5). In this restricted basis, the rotations 
properties (under <?_) are expressed by the fcs*fcs matrix Gt of coefficients 
j^f given by 119H (we define ks = l + 2[fc/10]). Thus, invariant functions 
are the eigenvector(s) of the matrix G&, corresponding to the eigenvalue 
1, when they exist. We note such an eigenvector by its fcs components 
(/m 2 )- (Note that, when the multiplicity of the eigenvalue 1 is larger than 
1, they can be several such vectors; this occurs for instance for the value 
k — 60, see the table.) Finally, 

the eigenfunctions of the dodecahedral space corresponding to the 
eigenvalue — k (k + 2) (k even) are all combinations of the functions 

J2 M 2 f M 2 1~k:M!,M 2 , 

where (/m 2 ) is an eigenvector of the matrix G^ with eigenvalue 1, 
for all the entire values of Mi between and k . 
For a value of k, the calculation of the (fcs) 2 coefficients of Gj, is 
immediate. The search for eigenvectors also runs very easily (on MAPLE 
for instance). Results confirms the values of the eigenvalues and their 
multiplicities given by |U. The table gives the numerical values of the 
first eigenvectors of the dodecahedron. MAPLE codes to calculate the 
modes may be obtained by request at marclr@cea.fr. 
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Table 1: For each value of k (from 12 to 44; eigenmodes exist only for 
the values given), the table give the fcs components (fcs = 1 + 2[fc/10]) of 
the vector /m 2 - The eigenmodes of D, corresponding to = — k (k + 2), 
are given by all combinations YIm f M 2 ^k-,M 1 ,M 2 ^ w here M2 varies from 
—k/2 and k/2 and verifies QP , and Mi takes all the entire values between 
-k/2 and k/2. 

k=12: .70352647068144845281, -.10050378152592120757 i, .70352647068144845294 

k= 20: .70702906968084266661, .010397486318835921579 i, 
-.0018904520579701675214, .010397486318835921658 i, .70702906968084267112 

k=24: .70697192267503596680, -.012403016187281332779 i, 
-.85867035142716919288 e-2, -.012403016187281332763 i, .70697192267503596704 

k= 30: .70709634963126361072+. 28263381117495907229 e-2 i, 
.10352886856225736179 e-4-.25900965187958376464 e-2 i, 
.23546331989053068040 e-3+.94117153238415742893 e-6 i, 0, 
-.23546331989053075280 e-3-.94117153238413699545 e-6 i, 
-.10352886856225582837 e-4+. 25900965187958382323 e-2 i, 
-.70709634963126377938-.28263381117495713938 e-2 i 

k= 32: .70708591611133287319, .54182828820791792693 e-2 i, 
.35336627491820729542 e-3, .22172001563495361997 e-3 i, 
.35336627491820731478 e-3, .54182828820791791292 e-2 i, 
.70708591611133284832 

k=36: .70708473425088583343, -.55056653389946616956 e-2 i, 
-.78652361985638027430 e-3, .70373165987149811020 e-3 i, 
-.78652361985638024036 e-3, -.55056653389946617208 e-2 i, 
.70708473425088584255 

k= 40: .70710661000066274387 , .49002537075583660713 e-3 i, 
-.44188504840885072726 e-4, -.39518175060953038014 e-5 i, 
.71851227383553307597 e-6, -.39518175060954735563 e-5 i, 
-.44188504840885377434 e-4, .49002537075583232202 e-3 i, 
.70710661000067363741 

k= 42: .70710034184850288923, -.30174979026251895955 e-2 i, 
.29583312770835158828e-4 -.19722208513890132964 e-4 i, 
.19722208513890135205 e-4 i, -.29583312770835191330 e-4 
.30174979026251928239 e-2 i, -.70710034184850303513 

k= 44: .70709762382971759036, .35958850229318122280 e-2 i, 
.13964602030803151449 e-3, -.15649985034520771577 e-4 i, 
.21878574396085900270 e-4, -.15649985034520774732 e-4 i, 
.13964602030803151858 e-3, .35958850229318113024 e-2 i, 
.70709762382971762555 
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Table 1, continued: idem, for values of k from 48 to 56 



k= 48: .70709753381495795208 -.36096867863497570254 e-2 i, 
-.20595810492686306199 e-3, .61892333738632857047 e-4 i, 
.56297546508021977011 e-4 .61892333738632854996 e-4 i, 
-.20595810492686303733 e-3, -.36096867863497572060 e-2 i, 
.70709753381495797378 



k=50: .70710677432308054872 , -.98118423819072486562 e-4 i, 
.88473375505798110701 e-5 .79771076275668425524 e-6 i, 
-.72519160250664405750 e-7 .13686767861056995799 e-19 i, 
.72519160250661739763 e-7, -.79771076275733196796 e-6 i, 
-.88473375505808361189 e-5, .98118423819152347705 e-4 i, 
-.70710677432345667134 



k=52: .70710632382180878363, .80371462823607746218e-3 i, 
-.29191620936716760517e-4, .12714219501894409823e-5 i, 
-.45691726334935131233e-6, -.14757576207556568245e-6 i, 
-.45691726334935055995e-6, .12714219501895343589e-5 i, 
-.29191620936716944313e-4, .80371462823612059084e-3 i, 
.70710632382183246593 



k= 54: .70709995013538603217, -.31081316489467430301 e-2 i, 
-.47140520459254381578 e-5, -.33908093663673232699 e-5 i, 
-.15713506819751032060 e-5, 0, 

.15713506819751014572 e-5 , .33908093663673239606 e-5 i, 
.47140520459254356244 e-5 , .31081316489467431476 e-2 i, 
-.70709995013538374766 



k=56: .70710165828128236325, .26907309019741289751 c-2 i, 
.69176592337987433371 e-4, -.59344090988012596170 e-5 i, 
-.63998529496876216108 c-6, -.20222732831614883787 e-5 i, 
-.63998529496876549920 e-6, -.59344090988012761721 e-5 i, 
.69176592337987085779 e-4, .26907309019741259726 e-2 i, 
.70710165828128410534 
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Table 1, continued: idem, for values of k from 60 to 62 



k= 60: 6.9394960301822869967-.50386496302764001912 i, 
.29358129642425900098-.64161677487791106593 i, 
-.25742682420776540173e- 2-. 11683763224239293660 e-2 i, 
-.36432453284660713864e-4+.77853455569667120805 e-4 i, 
.11085720861501942161e-4+.49965917891592449037 e-5 i, 
.21499250718118950501e-5-.46857540636509975228 e-5 i, 
-.39964490355299009740e-5-. 18268003343521902658 e-5 i, 
.21499250718136508514e-5-.46857540636526609287 e-5 i, 
.11085720861498572522e-4+.49965917891578880069 e-5 i, 
-.36432453284628003387e-4+.77853455569638907558 e-4 i, 
-.25742682420783301674e- 2-. 11683763224246853776 e-2 i, 
.29358129642426822804-.64161677487790595997 i, 
6.9394960301835349645-.50386496298031253450 i 

.80848705655829822065e-2+.70706052583861537173 i, 
-.20832969516412067927e-3+.65282867307451706578 e-4 i, 
.23918191153754173532e-6-.99923368048573892481 e-6 i, 
.18133484953835426669e-6-.98554877105180951336 e-8 i, 
-.94672012551906538764e-9+. 11022887144687929246 e-7 i, 
-.26567291467002942075e-8+.49083380949540372946 e-9 i, 
.40224546881829871331e-9+.94397550786428917031 e-9 i, 
-.26567291466915270685e-8+.49083380949346250026 e-9 i, 
-.94672012552856441064e-9+. 11022887144643159845 e-7 i, 
.18133484954122299052e-6-.98554877083976930526 e-8 i, 
.23918191155535933834e-6-.99923368050776437049 e-6 i, 
-.20832969516570937280e-3+.65282867306944956583 e-4 i, 
.80848705598891417065e-2+.70706052584512272298 i 

k= 62: .70710674764980394726, -.21751107313340935062 e-3 i, 
.10807153948145142931 e-4, .18045880549024964468 e-6 i, 
.55301892005028276440 e-7, .11642503580014562160 e-7 i, 0, 
-.11642503580010046955 e-7 i, -.55301892005094596638 e-7, 
-.18045880548966284252 e-6 i, -.10807153948148682328 e-4, 
.21751107313369991548 e-3 i, -.70710674765011395485] 
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